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ABSTRACT 

We analyse the u-r color distribution of 24346 galaxies with M,. < -1 8 and z < 0.08, drawn from the Sloan 
Digital Sky Survey first data release, as a function of luminosity and environment. The color distribution is 
well fit with two Gaussian distributions, which we use to divide the sample into a blue and red population. 
At fixed luminosity, the mean color of the blue (red) distribution is nearly independent of environment, with a 
weakly significant (~ 3cr) detection of a trend for colors to become redder by 0. 1-0. 14 (0.03-0.06) mag with a 
factor ~ 100 increase in local density, as characterised by the surface density of galaxies within a ±1000 km s~' 
redshift slice. In contrast, at fixed luminosity the fraction of galaxies in the red distribution is a strong function 
of local density, increasing from ~ 1 0-30 per cent of the population in the lowest density environments, to ^ 70 
per cent at the highest densities. The strength of this trend is similar for both the brightest (-23 < M,. < -22) 
and faintest (-19 <Mr< -18) galaxies in our sample. The fraction of red galaxies within the virialised regions 
of clusters shows no significant dependence on velocity dispersion. Even at the lowest densities explored, a 
substantial population of red galaxies exists, which might be fossil groups. We propose that most star-forming 
galaxies today evolve at a rate that is determined primarily by their intrinsic properties, and independent of 
their environment. Any environmentally triggered transformations from blue to red colors must either occur 
on a short timescale, or preferentially at high redshift, to preserve the simple Gaussian nature of the color 
distribution. The mechanism must be effective for both bright and faint galaxies. 
Subject headings: galaxies: clusters — galaxies: evolution 



L INTRODUCTION 

The local galaxy population is known to consist broadly 
of two types, identifiable for example by their morphologi- 
cal properties, and are termed late-type (spiral) and early-type 
(E/SO) galaxies. Recently, this division has been quantified 
in large da tasets, in the related quantities of broadband color 
JStrateva et al. 2001; Blanton et al. 2003a; Kauffmann et al. 
200 31) and star form ation rate (SFR, .Brinchmann et al...2003.; 



Balogh et al.l2004l hereafter Paper II). In particular, the color 



distribution at fixed luminosity is surprisingly well modelled 
by only two Gaussian distributions (Baldrv et al. 2004, here- 
after Paper I). The mean and variance of these two distribu- 
tions are strong functions of lumin osity, or stellar mass (see 
also Bernardi et al. 2003a; Blanto n et alJl200 3a: Hogg et aP 
12004); a similar trend is seen out to z ~ 1 JBell et al.i2004j) . 

Dressier ( 1980) was first to show that the fraction of late- 
type galaxies depends strongly on local galaxy environment, 
and rel ated trends in SFR have also been observed (e.g. Pa- 
per II, iMiller et aP 120031) . If these trends are due to an 
environmentally-induced transformation, we can look for 
clues to its nature in the change of the properties of galaxies 
within each type as a function of environment. For example, 
if the transformation is due to a decreasing SFR in late-type 
(blue) galaxies due to interactions with neighbouring galaxies, 
then their color distribution should be skewed toward redder 
colors in denser environments. 

Many studies of the fundamental plane have shown that the 
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properties of ear ly-type galaxies are nearly independent o f 
environment (e.g. lDressler et alJl987tlBernardi et alJ2003bl) . 
The properties of the blue population as a function of envi- 
ronment have been studied in less detail; early suggestions 
that their colors do not depend o n environment were seen in 
the small sample of lLarson et alJ ( 119 80). This result has been 
supported by analysis of the Tully-Fisher relation in clusters, 
groups and the field, using samples of several hundred galax- 
ies. In particular, the scatter in this relation is related to galaxy 
colo r, and it has been shown to be insensitiv e to environment 
re.g. lGinricinetalJl9 86: Bivi ano et alJl990l) . 

Using galaxy colors to describe the galaxy population, 
rather than morphological types, has the advantage that they 
are easily quantifiable, the measurements are robustly repro- 
ducible, and models exist which allow us to directly relate 
the m to star formation histor ies with minimal assumptions 
(e.g. iBruzual & CharloBl2003h . In particular, in Paper I we 
have fit the color distributions of g alaxies sele cted from the 
Sloan Digital Sky Survey (SDSS, Yorketal. 2000) with a 
two-component Gaussian model, as a function of luminosity. 
This technique has an advantage relative t o using strict color 
cuts to define separate populations (e.g. iButcher & Oemlen 
1984; Kodama & Bower 2001), as it accomodates photomet- 
ric and astrophysical uncertainties, such as aperture effects 
( B ernardi et al. 2003c). dust reddening and stochastic varia- 
tions in star-formation history. Furthermore, no a priori as- 
sumption about how to classify the population is required. 

In Paper I we used this novel procedure to show that the 
luminosity function of red galaxies is consistent with a model 
in which they are built up from mergers of blue galaxies. In 
this Letter, we develop this investigation by exploring how the 
model fits depend on local environment. 

2. data 

To allow us to reliably measure local galaxy densities, we 
select a sample of galaxies from the first data release of the 
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Fig. 1. — The solid points in each panel show the galaxy color distribution for the indicated 1 magnitude range of luminosity (right axis) and the range of local 
projected density, in Mpc"^, shown on the top axis. I-ct error bars are given by \/W+2, where N is the number of galaxies in each bin. The solid line is a double 
Gaussian model, with the dispersion of each distribution a function of luminosity only. The reduced x^ value of the fit is shown in each panel. 



SDSS f Abazaiia n et alJl20 03\ in the redshift range z < 0.08, 
with luminosities Mr < -18, assuming a ACDM cosmology 
with n,„ = 0.3, A = 0.7 and Ho = 70kms"' Mpc"'. Galaxy 
colors are measured from model magnitudes (IStoughton et al.l 
12902), corrected for Galactic extinction using the dust maps 
of Schlegel et al. ( 1998) and k-corrected to z = using the 
iBlanto n et al. (2003b) model. Projected local densities, E5, 
are computed from the distance to the fifth nearest neigh- 
bour that is brighter than M,- = -20 and within a redshift slice 
±1000kms~' of each galaxy, as described in Paper II. Apply- 
ing this bright limit to the density calculation gives us a uni- 
form density estimate that is applicable to our magnitude lim- 
ited sample over the full redshift range. To ensure robust mea- 
surements, we only consider galaxies sufficiently far from the 
survey boundary that the density estimate is unbiased, which 
limits our sample to 24346 galaxies. As a measurement of the 



environment on large scales, we use the C4 catalogue, based 
on the SDSS (Gomez et al. 2003, Miller et al. in prep), to 
identify galaxies that lie in groups or clusters. Cluster mem- 
bers are defined to be those within the virial radius and within 
1000 kms"' of the average cluster redshift. 

3. RESULTS 

FigureQ]shows the u-r color distribution of galaxies in bins 
of local density and luminosity. The population is divided 
into five density bins with the three middle bins having an 
equal number of galaxies and the least and most dense having 
half as many, to sample the extremes of the distribution. The 
median density of galaxies in the sample is E5 r^ 0.8 Mpc~^; 
thus, the lowest density bin (median density ~ 0.1 Mpc"^) is 
underdense by a factor ^ 8, while the densest bin corresponds 
to the typical density found in cluster cores. Following Pa- 
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per I, we model the distributions with two Gaussians using a 
Levenberg-Marquardt algorithm. The mean and amplitude of 
each distribution are varied as a function of luminosity and lo- 
cal density, while the dispersions are constrained to be a func- 
tion of luminosity only. The optimal parameters are found by 
minimizing the x^ fit to the data in Figure [0 these best-fit 
models are shown as the solid lines. 

The most striking result of Figure [2 is that the double- 
Gaussian model provides an acceptably good fit to all the 
distributions, with reduced x' ^ ^- Specifically, there is no 
convincing evidence that the shape of the blue distribution is 
significantly distorted in dense environments as might be ex- 
pected if these galaxies were slowly being transformed into 
the late-type (red) population. 

The dominant change in the galaxy population as a function 
of environment is in the relative number of galaxies in each 
distribution, as seen in earlier morphological (e.g. Dressier 
fl980i) and emission-line (e.g. Paper II) studies. To show 
this explicitly, we divide the sample into bins of luminosity 
and environment (either local density or group velocity dis- 
persion) and compute the relative abundance of red galaxies 
within each bin; the results are shown in Figure|2] There is a 
strong and continuous dependence on local density, with the 
fraction of galaxies in the red distribution at fixed luminosity 
increasing from 10-30 per cent of the population at the lowest 
densities, to ^ 70 per cent of the population in the highest den- 
sity environments. This is stronger than the dependence on lu- 
minosity at fixed density; in particular, the trend with density 
is of a similar magnitude at all luminosities. The fraction of 
red-distribution galaxies (at fixed luminosity) within the virial 
radius of clusters is independent of velocity dispersion, within 
the fairly large uncertainties. This implies that the population 
differences are primarily related to local galaxy density and 
not cluster mass or dynamics, in agreement with other work 



jDressleJlQSOllFairlev et al.l2002tlDe"Propris et a l.'200?). 

In contrast to the strong trend in the fraction of red galax- 
ies with E5, the mean color of each distribution depends only 
weakly on environment, as shown in Figure |3l Over a factor 
^ 100 in density, the mean color of the blue (red) popula- 
tion changes by only 0.1-0.14 (0.03-0.06) mag, depending 
on luminosity; the strongest trend is in the faint, blue galaxy 
distribution. These trends are weak, relative to the luminosity 
dependence: for example, the mean color of the blue popula- 
tion is ~ 0.7 mag redder in the brightest galaxies, compared 
with the faintest galaxies in the same environment. Similarly, 
we see little evidence that the dispersion depends on environ- 
ment: the reduced x^ of the fits do not improve significantly if 
the dispersions are allowed to vary, and there is no significant 
trend in the best-fit dispersions with environment. 

4. DISCUSSION AND CONCLUSIONS 

It has been known for a long time that the galaxy color dis- 
tribution (and correla ted quantities) depends on b oth environ- 
ment (e .g. Melnick & Sargeni 19771: Butcher & O emleJ 19781: 
[Dressier 1980; Mathieu & Soinrad 1981) and luminosity (e.g. 
|de Vaucouleurs 1961; Bower etal. 1992). We have shown 
here that the two effects are separable if the full color distri- 
bution is considered as the sum of two populations: while the 
mean and dispersion in the color of each type depend strongly 
on luminosity (Paper I), they are weak functions of environ- 
ment. For the red galaxy distribution, this is well known (e.g. 
ISandage & Visv anathan 1978). and is consistent with the idea 
that these are old galaxies, for which col or is only weakly 
sensitive to the present age jAndreonl2003ft . 
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Fig. 2. — The fraction of galaxies in the red distribution, in bins of local 
density (left panel) or cluster velocity dispersion (right panel), based on the 
Gaussian fitting of Fig. Q Five lines are shown, corresponding to five differ- 
ent luminosity ranges, as indicated. Points are shown with 1 - ct error bars 
at the median density or velocity dispersion of each bin. Each point repre- 
sents the fraction of red galaxies only among galaxies with that luminosity 
and environment and, therefore, they do not need to add to 100 per cent. 
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Fig. 3. — The mean color of the red distribution (left panel) and the blue 
distribution (right panel) as a function of local density, based on the Gaus- 
sian fitting of Fig.Q Five lines are shown, coiresponding to five different 
luminosity ranges, as indicated. Points with 1 — o" error bars are shown at the 
median density of each bin. 



However, it is a surprising result that the colors of blue 
galaxies, which are still actively growing and evolving, show 
such little dependence on environment, a result that was indi- 
rectly revealedthroughearlie£analy sis of the Tully-Fisher re- 
lation (e.g. Giuricin et al. 1986). For example, if the increased 
abundance of red galaxies is due to interactions (e.g. merg- 
ers and harassment) that increase in frequency monotonically 
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and smoothly with local galaxy density, we would expect the 
blue distribution to become increasingly non-Gaussian with 
density, and to gradually blend into the red distribution. It 
seems unlikely that the small change of '^ 0. 1 magnitudes in 
the mean color of the blue distribution as a function of en- 
vironment can be related to the much larger change in the 
abundance of this population, relative to the red population. 
Instead, it may indicate a small difference in the recent star 
formation history, or the distribution of intrinsic properties 
(e.g. dynamical mass or stellar velocity dispersion) of blue 
galaxies at fixed luminosity in different environments. 

We propose that characteristic properties (e.g. metallicity, 
dust content, and past-averaged SFR) of the late-type galaxies 
are determined primarily by their luminosity (likely through 
its relation to mass or other fundamental, intrinsic quantities) 
and that only interactions of a certain level trigger a trans- 
formation from late to early type. This transformation must 
be either sufficiently rapid, or sufficiently rare, to keep the 
overall color distribution unchanged. Furthermore, the mech- 
anism responsible for this transition must be effective for both 
bright and faint galaxies, as the trends with local density are 
of a similar magnitude for both populations. 

To quantify how rapid this transformation may be, we use 
the Bruzual & Chariot ( 2003 ) stellar population models to fol- 
low the color evolution from blue to red, assuming a range 
of exponential transformation timescales r. We start with a 
solar-metallicity model of a 9 Gyr old galaxy with a Salpeter 
( il955.) initial mass function, constant SFR, and a two com- 
ponent dust model with an effective extinction of 1 mag, of 
which 0.3 mag arises from the ambient interstellar medium 
|Chai-lot & Fall (2000). Brinchmann et al. (2003) have shown 
that this star formation history is a good description of the 
local, star forming population. We then assume the red popu- 
lation has been built up by a steady truncation of star forma- 
tion in blue galaxies, over a Hubble time. If the truncation 
is rapid (t < 0.5 Gyr), the color reddens to {u-r)o ~ 2.3 in 
only 0.75 Gyr. This change is so rapid that we expect only 
~ 1 per cent of the total galaxy population to be observed in 
the mid-transition phase, defined as the ~ 0.2 magnitude dip 
observed between the two Gaussian distributions. From Fig- 
ure Q it is evident that this would amount to an increase of 
only ^10 per cent in these intermediate-color bins, compara- 
ble to the observational uncertainties. Therefore, the observed 
color distribution is not significantly altered by the presence of 
this transforming population, and the simple Gaussian model 
remains a good fit to the data. 

On the other hand, if the SFR in a blue galaxy decays with 
with an exponential timescale of 2 Gyr, it takes ^ 4 Gyr for 



the colors to become as red as the early-type population. In 
this case, the number of galaxies with intermediate colors at 
the present day would increase by a factor ^ 2, and would 
distort the observed distribution in a way that is inconsistent 
with the observations. However, we caution that these results 
are sensitive to the assumption that the transformations occur 
uniformly in time. If, instead, transitions were more common 
in the past, the data can accomodate a slower rate of SFR 
decay, as expected if galaxies are stripped of their hot gas in 
dense environments (Larson et al. 1980; Balogh et al. 2000). 

We therefore conclude that short-timescale transformations 
could play a role at all densities and luminosities, without 
disrupting the Gaussian model fit. The small fraction of 
galaxies predicted to be in the transition phase is comparable 
to the fraction of spectroscopically identified post-starburst 
galaxies (Goto et al. 2003a; Ouintero et al. 2004), and anemic 
(passively evolving) spirals (van den Bergh 1976; Got o et alJ 
2003b) which may be the signature of such transformations 
(e.g. Dressier & Gunn 1992). 

It is interesting to compare the se resul ts with 
the morphology-density relation (iDressled Il980t 
Postman & Geller 1984!). In contrast with the bimodality of 
the color distribution, there are at least three morphological 
types that have different dependences on environment; 
ellipticals, SOs, and spirals. One possibility is that the 
transforming galaxies correspond to the SO population, since 
morphology may change on a longer tim escale than color i f 
the change is due to a decline in SFR (iBekki et al."2002!). 
This interpretation is in qualitative agreement with the dearth 
of such galaxies at higher redshift (Dressier e t al.l ll997i) . 
More detailed analysis will be the subject of future work, as 
morphologies for the SPSS galaxies become available (e.g. 
iKellv & McKavl2n04h . 

Finally, we note that, at all magnitudes, a population of red 
galaxies exists even in low-density environments; this is con- 
sistent with Paper II, where we found a population of galaxies 
without significant Ha emission in all environments. There- 
fore, either some fraction of the red population must arise in- 
dependently of environment (e.g., by consumption of the in- 
ternal gas supply), or these are fossil groups which result from 
the complete merging of bright galaxies (Ponman et.al,, 192a 
iMulchaev & Zabludoffll999HRomer et alJ2000h . 
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